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Abstract 


This paper presents a stationary and dynamic study of the advantages of using a regulating valve for the cathode outlet flow in combination with 
the compressor motor voltage as manipulated variables in a fuel cell system. At a given load current, the cathode input and output flowrate determine 
the cathode pressure and stoichiometry, and consequently determine the oxygen partial pressure, the generated voltage and the compressor power 
consumption. In order to maintain a high efficiency during operation, the cathode output regulating valve has to be adjusted to the operating 
conditions, specially marked by the current drawn from the stack. Besides, the appropriate valve manipulation produces an improvement in the 
transient response of the system. The influence of this input variable is exploited by implementing a predictive control strategy based on dynamic 
matrix control (DMC), using the compressor voltage and the cathode output regulating valve as manipulated variables. The objectives of this 
control strategy are to regulate both the fuel cell voltage and oxygen excess ratio in the cathode, and thus, to improve the system performance. All 


the simulation results have been obtained using the MATLAB-Simulink environment. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Polymer electrolyte membrane (PEM) fuel cell systems are 
efficient devices which allow the transformation of chemical 
energy stored in hydrogen to electric energy. In order to obtain 
this transformation efficiently the global system, which includes 
several subsystems, must be considered. The air supply is one of 
these subsystems, which has a great influence in the system effi- 
ciency. At the same time, one of the most important challenges 
in fuel cell control is to assure sufficient amount of oxygen in 
the cathode when current is drawn abruptly from the fuel cell 
stack. In this study, the air supply subsystem is composed by an 
air compressor, while the fuel supply subsystem relies on a pres- 
surized hydrogen tank. The hydrogen inlet flow rate is regulated 
by an independent control loop to maintain the working pressure 
in the anode close to the pressure in the cathode. A schematic 
diagram of the system is showed in Fig. 1. 
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The PEM fuel cell system without a proper controller will 
not be able to withstand the load fluctuations [1]. When an elec- 
tric load is connected to the fuel cell, the control system must 
maintain the optimal temperature, the membrane hydration and 
the partial pressure of the gases at both sides of the membrane 
in order to avoid voltage degradation and fuel cell life shorten- 
ing [2]. In particular, the air supply results critical in the system 
performance because the oxygen reacts instantly as current is 
drawn, whereas the oxygen supply is limited by the dynamics 
of the inlet manifold and the air compressor [3]. The operating 
air pressure and the air stoichiometric ratio provided to the stack 
by the air supply subsystem, control the oxygen partial pressure 
at the cathode catalyst layer. This partial pressure has major 
influence in the cathode polarization and thus, in the conversion 
efficiency [4]. 

In several publications, the control of air supply has been 
approached. In Ref. [5], it is demonstrated the convenience of 
regulating the oxygen excess ratio in the cathode to maximize 
the system net power. The oxygen excess ratio or stoichiomet- 
ric ratio is defined as the ratio of inlet oxygen flow to reacted 
oxygen flow in the cathode. In Refs. [2,5,6,7], the control of 
oxygen excess ratio is approached through the manipulation 
of the compressor motor voltage. In Refs. [2,5], feedforward 
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Nomenclature 


A fuel cell active area (cm?) 
At cathode output valve area (cm?) 
Cp specific heat capacity of air (1004Jkg~! K7!) 
E single cell open circuit voltage (V) 
f output free response 

F Faraday number (96,485 C) 

gi dynamic matrix coefficients 

G controller dynamic matrix 

AH higher heating value of hydrogen (286 kJ) 


i current density (A cm7?) 

I stack current (A) 

J cost function 

m control horizon 

Mo, oxygen molar mass (32 x 107°? kg mol” !) 
n number of cells in the stack 

p prediction horizon 


Po, PH, Oxygen and hydrogen partial pressure (Pa) 
P power (W) 

Q,R weight matrices 

t time (s) 

T temperature (K) 

u control variable 

U future control vector 
V voltage (V) 

w output setpoint 

W flow rate (kg s7!) 

ĵ predicted output 


Subscripts and superscripts 


atm atmospheric 

aux auxiliary 

CM compressor motor 
CP compressor 

FC fuel cell 

i,j index 

rct reacted 

sys system 

T transposed 


Greek letters 


y ratio of the specific heats of air (1.4) 
n efficiency 
ho» oxygen excess ratio 


control is employed, meanwhile in Refs. [6,7], model based 
predictive control (MPC) is employed. These works indicate 
that there is a severe conflict between oxygen excess control 
and net power dynamic response if no extra power source 
is employed. This limitation arises from the fact that all the 
power required by the air supply compressor comes from 
the fuel cell stack. In Ref. [6], supercapacitors are included 
as an ancillary power source in order to support the power 
peaks. 


The voltage applied to the compressor motor is a suitable 
manipulated variable in order to control the oxygen excess ratio 
in the cathode, as it is showed in Refs. [2,5,6,7]. In this paper, we 
propose to use a regulating valve for the cathode outlet flow in 
combination with the compressor motor voltage as manipulated 
variables to control both the oxygen excess ratio and the stack 
voltage. In Refs. [2,5,6,7], the stack output voltage is let unregu- 
lated. However, depending on the application, the regulation of 
the stack output voltage would be required. 

The paper is organized as follows: in Section 2, a description 
of the fuel cell system model employed is done. Stationary and 
transient analyses of system performance are done in Sections 3 
and 4, respectively. In Section 5, a control strategy based on pre- 
dictive control is proposed. The results are presented in Section 
6 and, finally, conclusions are stated in Section 7. 


2. Pem fuel cell system model description 


This work utilizes the PEM fuel cell model proposed by 
Pukrushpan et al. [2]. This nonlinear model is used to represent 
the transient behaviour of the plant in the simulation analysis. 
The control oriented model for automation applications includes 
the transient phenomenon of the compressor, the manifold filling 
dynamics (both anode and cathode), reactant partial pressures, 
and membrane humidity. The stack voltage predicted by the 
model depends on load current, partial pressure of hydrogen 
and oxygen, fuel cell temperature and the contents of water in 
the membrane. Spatial variations are not included and constant 
properties are assumed in all volumes. Only temporary varia- 
tions are present. The model also assumes that the inlet reactant 
flows in the cathode and in the anode can be humidified, heated 
and cooled rapidly. With respect to the considered dynamics, 
the model neglects the fast dynamics of the electrochemical 
reactions (time constant of 107! s). Temperature is treated as a 
constant parameter because its slow behaviour (time constant of 
10? s), allowing to be regulated by its own (slower) controller. 
The model represents a 75 kW fuel cell system with 381 cells. 
For more model details see Ref. [2]. 


3. Stationary analysis 


As it has been mentioned, the air supply subsystem has a 
strong influence in fuel cell performance. On the one hand, 
an insufficient air supply may cause oxygen starvation in the 
cathode, which causes voltage reduction and membrane life 
shortening. On the other hand, the compressor operation implies 
a power consumption that diminishes the system efficiency. The 
system efficiency is defined as [8]: 

Vre (Pre — Pax) 
Nsys = (1) 
1.482 Pec 
where the constant 1.482 results of using the higher heating value 
of hydrogen (AH = 286 kJ mol~ t ), Vec the fuel cell voltage, Prc 
the fuel cell generated power and Pay x is the power consumed by 
the auxiliary equipment. The compressor power represents the 
major consumption in the auxiliary subsystem and is the only 
parasitic power considered in this work. 
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Fig. 1. PEM fuel cell system. 


The operation at a higher pressure increments the generated 
voltage as a result of the increment in the cathode oxygen partial 
pressure and anode hydrogen partial pressure: 


E = 1.229 — 0.85 x 1073(Tpc — 298.15) 
+ 4.3085 x 107% Trc [In(pu,) + 0.5 In(po,)] (2) 


where E is the generated voltage in a single cell in open circuit 
in volts, Trc the fuel cell temperature in Kelvin and py, and po, 
are the cathode partial pressure of hydrogen and oxygen, respec- 
tively, in atm. Nevertheless, an increment in the operating pres- 
sure implies a higher consumption of power in the air compres- 
sor, contributing to a reduction in the system efficiency. In fact, 
an increment in the cathode pressure produces an increment in 
the supply manifold pressure and thus, an increment in the pres- 
sure ratio across the compressor and in the compressor power 
consumption. The power consumed by the air compressor is: 


y-1 


CpT; Eg 
Poy = 2 (=) 1| Wo 3) 
NCP Patm 


where Wcp is the compressor air flow rate, Pem the compressor 
power, Tatm the inlet air temperature in the compressor, ncp the 
compressor efficiency, Psm the supply manifold pressure, Cp 
the specific heat capacity of air, equal to 1004Jkg~! K7! and 
y is the ratio of the specific heats of air equal to 1.4. 

The compressor motor voltage as a control input allows 
regulating the oxygen partial pressure in the cathode. Aug- 
menting the compressor voltage the oxygen partial pressure 
increases. However, the compressor power consumption also 
increases. Thus, it is important to employ a complementary way 
to increment the oxygen partial pressure. A diminution of the 
area of the valve that closes the cathode air flow contributes to 
increase the cathode pressure and at the same time, contributes 
to decrease the air flow, the stoichiometry, and the oxygen 
concentration. The consumption of the compressor has two 
opposite trends: the trend to increase due to the pressure rise and 


the trend to decrease due to the flow reduction. When all these 
effects are taken into account, there is a positive balance in the 
total efficiency diminishing the valve area as we can observe 
analyzing the polarization (Fig. 2), efficiency (Fig. 3), power 
compressor consumption (Fig. 4) and oxygen partial pressure 
curves (Fig. 5) plotted for two different valve areas and a certain 
compressor voltage (Vem = 140 V). It is important to note that 
the increase in the efficiency is not for all current densities. 
When the current density is high the flow and concentration 
reduction have a greater influence than the pressure increase and 
the result is a decrease in the oxygen partial pressure. In effect, 
the oxygen consumption is higher at higher current densities: 


Mond, 
— ane 


4F g 


Wo, mS 


where Wo, rct is the rate of oxygen reacted, Mo, the molar 
mass of oxygen equal to 32 x 107°? kgmol™!, n is the cell’s 
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Fig. 2. Polarization curves with Vem = 140 V and different valve area. 
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Fig. 3. Efficiency curves with Vom = 140 V and different valve area. 
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Fig. 4. Compressor power consumption curves with Vem = 140 V and different 
valve area. 
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Fig. 5. Oxygen partial pressure curves with Vom =140 V and different valve 
area. 
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Fig. 6. Efficiency curves with Vem = 100 V and different valve area. 


number, A the active cell area equal to 381 cm?, F Faraday 
number equal to 96,485 C and i is the current density. Hence, 
the oxygen molar fraction in the cathode will be lower because 
of the oxygen starvation. 

Depending on the operating pressure, the increment in effi- 
ciency happens along different current density ranges. With 
low operating pressures (lower compressor voltage, Vom) the 
system efficiency increment due to the valve closure occurs 
only at low current densities, whereas with higher operating 
pressures (higher compressor voltage) the increment in sys- 
tem efficiency occurs along greater current densities ranges. 
This can be seen in Figs. 6—8. For Vom = 100 V, the efficiency 
increment occurs only for current densities below 0.34 A cm~? 
(Fig. 6), for Vom = 140 V, the increment occurs for current den- 
sities below 0.67 A cm~? (Fig. 3), while for Vem = 180 V, the 
increment occurs for al the current densities analysed (Fig. 7). 
As can be seen in Fig. 8 for Vem = 180 V, the oxygen partial 
pressure rise, as a result of diminishing the valve area, occurs at 
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Fig. 7. Efficiency curves with Vom = 180 V and different valve area. 
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Fig. 8. Oxygen partial pressure curves with Vem = 180 V and different valve 
area. 


all current densities. In Fig. 9, it can be seen how the efficiency 
changes with A; for values between 20 and 40 cm?. 

As a conclusion of this stationary analysis, it can be stated that 
the cathode output valve area as well as the compressor motor 
voltage have to be adjusted in order to have high efficiency. 


4. Transient analysis 


Besides the possible performance improvement observed 
in the stationary analysis there is also an improvement in 
the transient behaviour using the cathode air flow valve as a 
manipulated variable. A preliminary transient analysis is made 
employing a DMC control strategy (the details of this control 
strategy are explained in the next section). The controlled 
variable in this analysis is the stack voltage. A stack current 
disturbance from 200 to 210 A is applied to the fuel cell system 
in order to compare the disturbance rejection capability of the 


42 


40 


wo 
œ 


Efficiency [%] 
8 


wo 
A 


—¥— A =350m" 


32 n 
a A, =40cm 


pm ee een penn Sl 


| 
0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Current density [A cm] 


Fig. 9. Efficiency curves with Vom = 180 V and different valve area. 
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Fig. 10. Comparative between stack voltage responses. 


system. A comparison between the transient responses obtained 
using only the compressor voltage as a control variable and the 
one obtained by using the opening valve area in combination 
with the compressor voltage is showed in Fig. 10. The later 
shows a better behaviour with a reduction in the stack voltage 
time response from 0.6 to 0.4s. 


5. Control strategy 


The advantages of using the opening valve area together with 
the compressor voltage are exploited implementing a control 
strategy. This control strategy is based in dynamic matrix predic- 
tive control (DMC) and has two control objectives: (i) to regulate 
the oxygen excess ratio in the cathode, Ao, and (ii) the generated 
voltage, Vac. As we explained before, the control of Ao, is an 
indirect way to control the system efficiency [5]. The load cur- 
rent, Ipc, is considered as a disturbance to the fuel cell system. 
The compressor motor voltage by itself is not able to control 
both control objectives. Nevertheless, with the addition of the 
opening valve area control variable, A;, this problem is solved. 
Therefore, the manipulated variables are: uj = Vem and uz = At. 

The DMC uses the step response to model the process, taking 
into account only the p first samples until the response tends to 
a constant value, assuming therefore that the process is asymp- 
totically stable. Thus, the predicted output can be expressed as: 


$= Gu+ f (5) 
where 
81 0 0 
82 = Bl 0 
G = (6) 
&m S8m-1 °*"° §1 
Sp Sp-1 §p—m+1 
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is the dynamic matrix constructed from the coefficients g; 
obtained from the step response with prediction horizon p and 
control horizon m, u the future control vector and f is free 
response, that means, the response that not depends of future 
control movements [9]. 

The DMC controller objective is to minimize the difference 
between the consignees, w(t), and the process outputs in a least 
square sense with the possibility of including a penalty term on 
the control signal: 


Pp m 
J =X REH jD- we + AP +Y olut j- DP 
j=l j=1 
(7) 


where R and Q are diagonal weight matrixes. With the matrix R 
it is possible to compensate the different ranges of values in the 
process outputs, meanwhile the matrix Q allows to give different 
weights to the control signals. 

If there are no restrictions in the manipulated variables, the 
minimization of the cost function J can be realized making its 
derivative equal to zero, resulting: 


U =(G™RG+Q) | G™RW- f) (8) 


Such as all the predictive control strategies, only the first ele- 
ment in the control vector calculated is sent to the plant. In the 
next iteration, the sequence of control is calculated again using 
actualized information from the plant. The theoretical bases of 
the method are given in detail in Ref. [9]. 

There is a great interaction between the manipulated variables 
and the controlled variables, which makes difficult the realiza- 
tion of a decentralized control of the system [10]. The DMC 
resolves efficiently the interaction problem between manip- 
ulated variables and controlled variables. The multivariable 
controller implemented employs in his control algorithm an 
extended dynamic matrix that takes into account the interactions: 


G G 

G= 11 12 (9) 
G21 G22 

where each matrix Gj contains the coefficients of the ith step 


response corresponding to the jth input. The matrices R and Q 
have now the following form: 


te DY. (a 0 
Eo G) 


where R; and Q; are diagonal matrices of dimension p x p and 
m x m, respectively. Thus, a centralized multivariable controller 
it is proposed, which take into account the interactions between 
manipulated and controlled variables. 


6. Results 


The control horizon m, the prediction control p, and the 
matrices R and Q were adjusted, in the control algorithm, to 
achieve an adequate dynamic behaviour of the fuel cell system. 
The values of R and Q have a strong influence on the transient 
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Fig. 11. Fuel cell stack current. 


response obtained. This is especially true for the weight matrix 
Q, which reduces the control effort. The higher are the values 
in Q, the lower is the control effort, but the response time is 
greater. The values chosen of Q and R are: Qı =diag(1,1, ..., 
1), Q2 = diag(0.5,0.5, ..., 0.5), Ry =diag(5,5, ..., 5), Ro=(.,1, 
..., 1). The values of m and p are 100 and 15, respectively. 

The internal model utilized by the DMC controller has 
been obtained through the linearization of the nonlinear model 
described in Section 2 around the operating point corresponding 
Vom = 187.5 V, Ar = 20cm, Ig, = 190 A, which gives a net power 
of 40 kW. 

The simulation results of the controlled system with varia- 
tions in the load current are presented in Figs. 12-15. In Fig. 11, 
the perturbation profile is showed. Figs. 12 and 13 show the 
controlled variables (Ag, and Vgc), while Figs. 14 and 15 show 
the manipulated variables (Vem and A;). As can be observed in 
the figures, the implemented controller has a good disturbance 
rejection: the stack voltage response has a peak (3.5% of Vac 
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Fig. 12. Oxygen excess ratio. 
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Fig. 13. Fuel cell stack voltage. 
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Fig. 14. Compressor motor voltage. 


Fig. 15. Cathode air flow opening valve area. 
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Fig. 16. Comparative between stack voltage response with measured and non- 
measured disturbance. 


variation for 15% of Ipc variation) that disappears in less than 
0.5 s, and the oxygen excess ratio response presents a peak (30% 
of Ao, variation for the same Jpc variation) that vanishes in less 
than 1s. 

The controller performance can be improved if the distur- 
bance is measured and a step disturbance model is included in 
the control model used to predict the plant outputs. Comparisons 
between measured and non-measured disturbance approaches 
have been made. In Figs. 16 and 17, details of voltage stack 
and oxygen excess ratio responses are showed, respectively. The 
results obtained show a substantial improvement in the gener- 
ated voltage response. The oxygen excess ratio response could 
be changed by adjusting the values of the weight matrices. 

The control algorithm implements a compensation mech- 
anism to rectify the inevitable model errors and deal with 
non-measured disturbances. This compensation mechanism uti- 
lizes output values, and thus, assures zero error at steady-state. 


]-t---4--4 


p 
HI 
A 
[i 
1 
! 
L 
i 
[i 
k 
1 
E 
19 195 20 205 21 215 22 225 23 235 24 
Time (s) 


Fig. 17. Comparative between oxygen excess ratio response with measured and 
non-measured disturbance. 
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The simulation analysis shows a good performance in a wide 
operating range around the linearization point despite the inter- 
nal controller model is linear. 


7. Conclusions 


High efficiency level, long life cycle and good transient 
behaviour are fundamental issues for the success of fuel cell 
systems in energy and automotive market. Through a steady- 
state analysis it is showed that the system efficiency, in most of 
the current density range, can be improved by manipulating the 
cathode outlet air flow valve. The dynamic analysis performed 
also shows a transient response improvement with this addi- 
tional manipulated variable. Taking advantage of these facts, 
this article proposes a control strategy based on predictive con- 
trol (DMC) that uses the compressor motor voltage together 
with the cathode air flow valve area as manipulated variables. 
The controlled variables are the stack voltage and the oxygen 
excess ratio. To predict the future process response, the control 
strategy makes use of a process model that is easily obtainable 
through step response. The simulation analysis shows an appro- 
priate dynamic response, which can still be improved with the 
inclusion of the disturbance model. The control objectives have 
been accomplished with reduced control effort. This effort can 
be further reduced modifying the values of the matrix Q. This 
is particularly important because of practical limitations in the 
manipulated variables. 
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